ASYMPTOTIC STABILITY OF THE PERIODIC TODA LATTICE 
UNDER SHORT-RANGE PERTURBATIONS 



SPYRIDON KAMVISSIS AND GERALD TESCHL 



Abstract. We consider the question of the asymptotic stability of periodic 
(and slightly more generally of algebro-geometric finite-gap) solutions of the 
doubly infinite Toda lattice under a short-range perturbation. We prove that 
the perturbed lattice asymptotically approaches a modulated lattice instead. 

More precisely, let g be the genus of the hyperelliptic curve associated with 
the unperturbed solution. We show that, apart from the phenomenon of the 
solitons travelling on the quasi-periodic background, the n/t-pane contains 
g + 2 areas where the perturbed solution is close to a finite-gap solution on 
the same isospectral torus. In between there are g + 1 regions where the per- 
turbed solution is asymptotically close to a modulated lattice which undergoes 
a continuous phase transition (in the Jacobian variety) and which interpolates 
between these isospectral solutions. In the special case of the free lattice 
{g = 0) the isospectral torus consists of just one point and we recover the 
known result. 

Both the solutions in the isospectral torus and the phase transition are 
explicitly characterized in terms of Abelian integrals on the underlying hyper- 
elliptic curve. 

Our method relies on the equivalence of the inverse spectral problem to a 
vector Riemann— Hilbert problem defined on the hyperelliptic curve and gen- 
eralizes the so-called nonlinear stationary phase/steepest descent method for 
Riemann-Hilbert problem deformations to Riemann surfaces. 



1. Introduction 

A classical result going back to Zabusky and Kruskal [IS] states that a decaying 
(fast enough) perturbation of the constant solution of a soliton equation eventually 
splits into a number of " solitons" : localized travelling waves that preserve their 
shape and velocity after interaction. This is the motivation for the result presented 
here. Our aim is to investigate the case where the constant background solution 
is replaced by a periodic one. We provide the detailed analysis in the case of the 
Toda lattice though it is clear that our methods apply to other soliton equations 
as well. 

The classical result is that an arbitrary short-range solution will eventually 
split into a number of solitons plus a decaying radiation part. In the case of the 
Korteweg-de Vries equation this was first proven by Sabat |,36j and by Tanaka _3S] . 
Precise asymptotics for the radiation part were first formally derived by Zakharov 
and Manakov [4? and by Ablowitz and Segur T, fST" with further extensions by 
Buslaev and Sukhanov 5^. A detailed rigorous justification not requiring any a 
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Figure 1. Numerically computed solution of the Toda lattice, 
with initial condition a period two solution perturbed at one point 
in the middle. 

priori information on the asymptotic form of the solution was given by Deift and 
Zhou [6] inspired by earlier work of Manakov ^ and Its [18] (see also [19], [20] . 
[21]). For further information on the history of this problem we refer to the survey 
by Deift, Its, and Zhou 

A naive guess would be that the perturbed periodic lattice approaches the un- 
perturbed one in the uniform norm. However, as pointed out in [24] this is wrong: 
In Figure [1] the two observed lines express the variables a{n,t) of the Toda lattice 
(see (jl.ip below) at a frozen time t. In areas where the lines seem to be continuous 
this is due to the fact that we have plotted a huge number of particles and also 
due to the 2-periodicity in space. So one can think of the two lines as the even- 
and odd-numbered particles of the lattice. We first note the single soliton which 
separates two regions of apparent periodicity on the left. Also, after the soliton, 
we observe three different areas with apparently periodic solutions of period two. 
Finally there are some transitional regions in between which interpolate between 
the different period two regions. It is the purpose of this paper to give a rigorous 
and complete mathematical explanation of this picture. This will be done by for- 
mulating the inverse spectral problem as a vector Riemann-Hilbert problem on the 
underlying hyperelliptic curve and extending the nonlinear steepest descent method 
to this new setting. 

Consider the doubly infinite Toda lattice in Flaschka's variables (see e.g. [14] . 

kn,t)^2{a{n,tf-a{n-l,tf), 
a(n, t) — a(n, t){b{n + l,t) — h{n, t)), 

(n, t) E Z X M., where the dot denotes differentiation with respect to time. 

In case of a constant background the long-time asymptotics were first computed 
by Novokshenov and Habibullin [33] and later made rigorous by Kamvissis [22] 
under the the additional assumption that no solitons are present. The case with 
solitons was recently solved by Kriiger and Teschl in [27] (for a review see [28]). 

Here we will consider a quasi-periodic algebro-geometric background solution 
{aq,bq), to be described in the next section, plus a short-range perturbation (a, &) 
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satisfying 

(1.2) ^|n|6(|a(n,i) — aq{n, t)\ + \b{n, t) — bq{n, t)\) < oo 

n 

for one (and hence for all, see [11]) t € M. The perturbed solution can be computed 
via the inverse scattering transform. The case where (a^, bq) is constant is classical 
(see again [M], [40] or [43]), while the more general case we want here was solved 
only recently in [11] (see also [31]). 

To fix our background solution, consider a hyperelliptic Riemann surface of genus 
g with real moduli Eq,Ei, ....,i?2g+i. Choose a Dirichlet divisor Vp^ and introduce 

(1.3) z(n, t) = 4„((X3+) - &p„{V^) - nA^_ (00+) + tH, - S^^ G 

where A^^^ fepo) ^'^ Abel's map (for divisors) and S^^^ , J7q are some constants defined 
in Section [2] Then our background solution is given in terms of Riemann theta 
functions (defined in (2.14)) by 

,,2 _ .2 ^{z{n + l,t))e{z{n-l,t)) 
"^^'"'^^ -'^ 6{z{n,t)Y 

(1.4) bqin,t) = b+--\og' " 



2dt \e{z{n - l,t)) 

where a, 6 e R are again some constants. 

Assume for simplicity that the Jacobi operator 

(1.5) H{t)f{n) = a(n, t)f{n + 1) + a{n - 1, t)f{n - 1) + b(n, t)f{n), f G i^Z), 

corresponding to the perturbed problem (jl.ip has no eigenvalues. In this paper we 
prove that for long times the perturbed Toda lattice is asymptotically close to the 
following limiting lattice defined by 



ai{j,t) 2 OUin,t)) e{z{n-l,t)+6{n,t)) 
^J-'aq{j,t)' e{z{n-l,t)) e{z{n,t)+S{n,t)) 



X 



(1-6) xexp(— / log(l-|i?nc.^^_ 

'C(n/t) 




5i{n,t)=^[ log(l - |i?nc,, 

lc{n/t) 



2ti\ 



where R is the associated reflection coefhcient, C« a canonical basis of holomorphic 
differentials, Woo+oo_ is an Abelian differential of the third kind defined in (j2.15p . 
C{n/t) — ■K~^{a{Hq) n (—00, Zj{n/t))) oriented such that the upper sheet is to the 
left, and Zj{n/t) is a special stationary phase point for the phase defined in the 
beginning of Section |4j From the products above, one easily recovers ai{n, t). More 
precisely, we have the following. 

Theorem 1.1. Let C be any (large) positive number and S be any (small) positive 
number. Let <E S be the 'resonance points' defined by S — {Eg : \R{Es)\ = 1}. 
(There are at most 2g + 2 such points, since they are always endpoints Ej of the 
bands that constitute the spectrum of the Jacobi operator.) Consider the region 
D = {{n,t) : If I < C} n {{n,t) : \zj{f) - E^l > 6}, where Zj(f ) is the special 
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stationary phase point for the phase defined in the beginning of Section Then 
one has 



(1.7) 



n 



1 



uniformly in D, as t ^ oo. 

The proof of this theorem wiU be given in Section 2] of this paper. 

Remark 1.2. (i) If eigenvalues are present we can apply appropriate Darboux 
transformations to add the effect of such eigenvalues ('[13] What we then see 
asymptotically is travelling solitons in a periodic background. Note that this will 
change the asymptotics on one side. In any case, our method works unaltered for 
such situations (cf. |12| ) as well. 

(a) Employing the very same methods of the paper it is very easy to show that in 
any region \j\ > C, one has 



(1.8) 



n 



a{i,t) 



1 



uniformly in t, as n ^ oo. 

(Hi) The effect of the resonances Eg is only felt locally (and to higher order in 1/t) 
in some small (decaying as t oo) region, where in fact \zj{j) — Es\ — > as 
t — >■ cxD. We expect a " collisionless shock" phenomenon to appear (|9], |23| ). We 
will study this in a future paper. 

(iv) For the proof of Theorem \l.l\ and Theorem \1.3\ it would suffice to assume ()1.2|) 
with \n\^ replaced by |np (or even \n\ plus the requirement that the associated re- 
flection coefficient is Holder continuous). Our stronger assumption is only required 
for the detailed decay estimates in Theorem \1.4\ below. 



By dividing in (|1.6p one recovers the a{n,t). It foUows from the main Theorem 
and the last remark above that 

(1.9) \a{n,t) - ai{n,t)\ ^ 

uniformly in D, as f — > oo. In other words, the perturbed Toda lattice is asymp- 
totically close to the limiting lattice above. 

A similar theorem can be proved for the velocities b{n,t). 



Theorem 1.3. In the region D 
of Theorem \1.1\ we also have 



{{n,t) : \f\<C}n{{n,t) : \z,m~Eg\>6}, 



Y,{biij,t)^b,{j,t)) 



(1.10) 

uniformly in D, as t ^ oo, where bi is given by 

oo 

E - b,{j,t)) =— / log(l - \R\^)no 

t^. ^TTl Jc(n/t) 



(1.11) 



ld_ 
2ds 



log 



0{z{n,s) +S{n,t)) 



and Qq is an Abelian differential of the second kind defined in (|2.16p . 
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The proof of this theorem wiU also be given in Section |4] of this paper. 

The next question we address here concerns the higher order asymptotics. Namely, 
what is the rate at which the perturbed lattice approaches the limiting lattice? Even 
more, what is the exact asymptotic formula? 

Theorem 1.4. Let Dj be the sector Dj — {{n, t), : Zj{n/t) G [E2j + e, £'2^+1 — e] 
for some £ > 0. Then one has 



(1.13) Y. Wi.0-'Ki,0) = y'^;7p-^;^2Re(«.>.i)iAi(.>.f))+O(r'') 



for any a < 1 uniformly in Dj, as i —> 00. Here 

(1.14) ^"(z,)/i=nuA^if^>o, 

1-^25+2(2^) 

1 /2 

(where (j){p^n/t) is the phase function defined in (|3.17p and R2g-\.2{^) square 
root of the underlying Riemann surface), 

/OO- 
^Ut{n.t),oCk{Zj), 

(1.15) Ai(n,t) = Woo_,o(^j) - X!'^'=^(^("'^))'^^><!(»^*)^o(oo+)Cfe(2;j)> 



with Cki{v{n,t)) some constants defined in (|5.14p . Wg.o an Abelian differential of 
the second kind with a second order pole at q (cf. Remark \5.3\) . 



/3 =^e'*''/'^""S(fl.(zj)))+arg(r(ii.))-2i.a(2:j)) ^ '/'"(^j) ^ ^-t4>{z,) ^ 

9{z{z,,n,t) + 6in,t)) ^^(^(^*, 0, 0)) 



6'(z(zj,0,0)) 0{z{z*,n,t)+6{n,t)) 

( ^ f ( 1-Ii?'^ 

(1.16) X exp — / log 



^27ri^(„/,/-Vl-|i?(-.)IV^'7' 
where r(z) is the gamma function, 

(1.17) ^=_j_iog(i_|i?(z,)n>o, 

and OLi^Zj) is a constant defined in (j4.22l) . 

The proof of this theorem will be given in Section [S] of this paper. 

Remark 1.5. (i) Combining our technique with the one from [7 enables one 

to find a complete asymptotic expansion. 
(ii) The necessary changes needed to include solitons are given in |29j (see also 
[13], [27], and [42];. 
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(iii) Finally, we note that the same proof works even if there are different spatial 
asymptotics as n ^ ±00 as long as they lie in the same isospectral class 

(cf m)- 

2. Algebro-geometric quasi-periodic finite-gap solutions 

As a preparation we need some facts on our background solution (og, hq) which 
we want to choose from the class of algebro-geometric quasi-periodic finite-gap 
solutions, that is the class of stationary solutions of the Toda hierarchy, [3], |16j . 
|40) . In particular, this class contains all periodic solutions. We will use the same 
notation as in |40| . where we also refer to for proofs. As a reference for Riemann 
surfaces in this context we recommend [15J. 

To set the stage let M be the Riemann surface associated with the following 
function 

25-1-1 

(2.1) Rl'gl^iz), R2g+2{z)=\{{z^Ej), Eo<E,<---<E2g+u 

g E N. M is a compact, hyperelliptic Riemann surface of genus g. We will choose 
^2g+2(^) the fixed branch 

2g+i 

(2.2) <+2(^) = - n V^^' 
where ^ is the standard root with branch cut along (— oo,0). 

1 /2 

A point on M is denoted hy p = (z, ±i?2g+2(^)) — {z , ±) , z E C, or p — (00, ±) = 
oo±, and the projection onto C U {00} by 7r(p) = z. The points {{Ej, 0),Q < j < 
2^ + 1} ^ M are called branch points and the sets 

(2.3) n± = {(z,±i?^/^2(z)) I z € C \ (j[E2j,E2,+i]} C M 

are called upper, lower sheet, respectively. 

Let {oj, bj}j^^ be loops on the surface M representing the canonical generators 
of the fundamental group 7ri(M[). We require aj to surround the points E2j^i, i?2j 
(thereby changing sheets twice) and bj to surround Eq, -B2J-1 counterclockwise on 
the upper sheet, with pairwise intersection indices given by 

(2.4) Gi o aj — bi o bj — 0, Oi o bj — 6ij, ^ l£ i, j l£ g- 

The corresponding canonical basis {Cj}j=i for the space of holomorphic differentials 
can be constructed by 

(2-5) C = EeW^i7^' 

where the constants c(.) are given by 

(2.6) c,(.).C- ^ = 2/ ^^eM. 
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The differentials fulfill 

(2.7) / Ck = Sj.k, / Ck = Tj^k, Tj,k^Tkj, l<j,k<g. 

J aj Jbj 

Now pick g numbers (the Dirichlet eigenvalues) 

(2.8) (A,)?=i-(A^„^,)f=i 

whose projections lie in the spectral gaps, that is, fij e [E2j~i, E2j]. Associated 
with these numbers is the divisor I?^ which is one at the points fij and zero else. 
Using this divisor we introduce 

z{p, n, t) = (p) - ttp^ (Vp) - nA^_ (oo+) + tUo - E^^^ G C^, 

(2.9) z(n, t) = z(oo+, n, t), 

where S^^ is the vector of Riemann constants 

(2.10) = i^^l^i^, Po = iEo,0), 

U_Q are the 6-periods of the Abelian differential Hq defined below, and A^^ (Sp^) is 
Abel's map (for divisors). The hat indicates that we regard it as a (single- valued) 
map from M (the fundamental polygon associated with M by cutting along the a 
and b cycles) to C^. We recall that the function 9{z{p,n,t)) has precisely g zeros 
fj.j{n, t) (with flj{0, 0) = fij), where 9{z) is the Riemann theta function of M. 
Then our background solution is given by 

^ e{z{n,t))^ 

^'■''^ '^^^^'^-'+2dM eizin-lt)) )- 

The constants a, b depend only on the Riemann surface (see [40l Section 9.2]). 
Introduce the time dependent Baker- Akhiezer function 

(2.12) ^,{p,n,t) ^ g(»,0,^) l^-^"'^ii exp (n [\oo^oo^+t Tl^o), 

6'(aP)0,o)) V j,^^ + j,^^ J 

where C(n, 0, t) is real-valued. 



2_ %(O,O))0(z(-l,O)) 



(2.13) C(n,0,t) 

a{z{n, t))6[z_(n — 1, t)) 

and the sign has to be chosen in accordance with aq(n, t). Here 

(2.14) e{z) = 2m ({m,z) + Ml^M^J ^ z € C", 

is the Riemann theta function associated with M, 

^n^K\ nf=i(i'-Aj), 

(2.15) Woo+oo_ = Yj2 "''^ 

is the Abelian differential of the third kind with poles at oo^ and oo_ and 

(2.16) llo = ^^i^rf., EA. = ^E^.' 
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is the Abelian differential of the second kind with second order poles at 00+ respec- 
tively oo_ (see [ini Sects. 13.1, 13.2]). All Abelian differentials are normalized to 
have vanishing aj periods. 

The Baker- Akhiezer function is a meromorphic function on M \ {oo±} with an 
essential singularity aX oo±. The two branches are denoted by 

(2.17) ijg^±{z,n,t) = tpg{p,n,t), p={z,±) 
and it satisfies 

(2.18) ^^^^{p^n,t) = P5,2(t)V'?(p,",i), 
where 

(2.19) Hq{t)f{n) = a,(n, t)/(n + 1) + a,(n - 1, t)f{n - 1) + 6,(71, t)f{n), 

(2.20) Pg,2(t)/(n) = a,(n, t)f{n + 1) - aq{n - 1, t)f{n - 1) 

are the operators from the Lax pair for the Toda lattice. 

It is well known that the spectrum of Hq{t) is time independent and consists of 
(7 -I- 1 bands 

9 

(2.21) a{Hq)^[j[E2,,E2,+i]. 

For further information and proofs we refer to [40, Chap. 9 and Sect. 13.2]. 
3. The Inverse scattering transform and the Riemann-Hilbert 

PROBLEM 

In this section our notation and results are taken from [10] and [11]. Let 
'09, ±(-2^7 n, t) be the branches of the Baker- Akhiezer function defined in the previous 
section. Let ip±{z,n,t) be the Jost functions for the perturbed problem 
(3.1) 

a{n, t)ip±{z, n + l,t) + a{n — 1, t)'ip±{z, n — l,t) + b{n, t)'il;±{z, n, t) = z'ip±{z, n, t) 
defined by the asymptotic normalization 

(3.2) lim w(z)T"(Vj±(z,ri,t)-V,.±(z,n,<)) = 0, 
where w{z) is the quasimomentum map 

(3.3) w{z)=exp{f ujoo+00^), p^{z,+). 

J Eq 

The asymptotics of the two projections of the Jost function are 



z 

ip±{z,n,t) = V'g,±(z,0,t) 



(3.4) X (l+(i?±(n,i)±^6,(j--«,t))- + 0(-)), 
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as z — 7- OO, where 



(3.5) 



j=n 3=n+l 

n-1 



A-{n,t)=l[ B.{n,t)= J2 Mj,t) - b{j,t)). 



j = -co 

One has the scattering relations 



(3.6) T{z)->p^{z,n,t) = ip±{z,n,t) + R±iz)'ip±{z,n,t), z e a{Hg), 

where T(z), R±{z) are the transmission respectively reflection coefficients. Here 

ijj±{z,n,t) is defined such that ■4'±{z,n,t) = \\mg:^Qtl)±{z -\- ie,n,t), z G <7{Hq). If 
we take the limit from the other side we have 'tp±{z^ n, t) = lime4.o '0±(-z — i£, n, t). 
The transmission T{z) and reflection R±{z) coefficients satisfy 



(3.7) T{z)R+{z)+T{z)R-{z) = Q, |T(z)|^ + |i?±(z)|^ = 1. 

In particular one reflection coefficient, say R{z) = i?+(z), suffices. 

We will define a Riemann-Hilbert problem on the Riemann surface M as follows: 

(3.«j m(p,n,^j_| (V,+ (z,n,t) T{z)i^.{z,n,t)) , p={z,-) ' 

Note that m{p,n,t) inherits the poles at /ij(0, 0) and the essential singularity at 
oo± from the Baker-Akhiezer function. 

We are interested in the jump condition of m{p, n, t) on E, the boundary of II-i- 
(oriented counterclockwise when viewed from top sheet n+). It consists of two 
copies S± of <j{Hq) which correspond to non-tangential limits from p = {z, +) with 
±Im(2;) > 0, respectively to non-tangential limits from p — {z,—) with =Flm(2;) > 0. 

To formulate our jump condition we use the following convention: When repre- 
senting functions on E, the lower subscript denotes the non-tangential limit from 
n_|_ or n_, respectively, 

(3.9) rn±(po) = lim m{p), po € E. 

Il±3p->-po 

Using the notation above implicitly assumes that these limits exist in the sense that 

m{p) extends to a continuous function on the boundary away from the band edges. 

Moreover, we will also use symmetries with respect to the the sheet exchange 
map 



(3.10) p* 
and complex conjugation 

(3.11) p={ 




for p = {z,±), 
for p = oo±, 



{z,±) forp=(^,±)^E, 
{z,t) {oip={z,±)eT,, 
oo± forp = oo±. 



In particular, we have p = p* for p G E. 

Note that we have m±(p) = m;p(p*) for rh{p) = m(p*) (since * reverses the 
orientation of E) and rh±{p) = m±{p*) for m{p) = m{p). 
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With this notation, using p.6p and p.7|) . we obtain 



(3.12) m+{p,n,t) = m^{p,n,t) (^'^(j]'' ^^^P^^ , 

where we have extended our definition of T to E such that it is equal to T{z) on 
S+ and equal to T(z) on E_. Similarly for R{z). In particular, the condition on 
S+ is just the complex conjugate of the one on E_ since we have R(j>*) = R{p) 
and m±{p* , n, t) = m±{p, n, t) for p e S. 

To remove the essential singularity at oo± and to get a meromorphic Riemann- 
Hilbert problem we set 

(3.13) m^ip,n,t) = mip,n,t)(^^^^P*';;^'^-' ■ 
Its divisor satisfies 

(3.14) (m?) > -I5^(„,4)., (m^) > -2)^(„,4), 
and the jump conditions become 

ni'^{p, n, t) = 'rrL^{p, n, t)j'^{p, n, t) 



r315l j'(nnt)-( 1-1^(^)1' ~RipMp,n,t)e-'^(PA 

(3.15) ^ l^i?(p)e(p,n,Oe**(p) 1 J 

9{zip,n,t))0iz{p*,O,O)) 



where 

(3.16) &{p,n,t) 



9{zip,0,0))e{z{p*,n,t)) 



and 



(3.17) <P{p,-)^2 no + 2- cjoo+oo_ei 



for p e E. Note 

Observe that 
and 



#^-e(,,n,.)e*^(^). 



m (p) — m?{p) 



1 

1 



TO {p*) — m [p) 

which follow directly from the definition (13. 13^ . They are related to the symmetries 

J\P)=1W) and J\p)=(^^ J) •^'(^*)"' (! J)- 

Now we come to the normalization condition at cx)+. To this end note 
(3.18) 

m{p,n,t)= (^A+{n,t){l-B+{n~l,t)^) -^-i-^(l + i3+(n, t)i)) + 0( J_), 

for p = {z,+) — > oo-|_, with A±{n,t) and B±{n,t) are defined in p.Sp . The formula 

near oo- follows by flipping the columns. Here we have used 

(3.19) 

r(z) = A_(n,t)A+(n,t)(l - B+{n,t) + h,M-h{n,t) + BAn.t) ^ 
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Using the properties of '>p{p,n,t) and '>jjq{p,n,t) one checks that its divisor satisfies 

(3.20) (toi) > -2?^(„,t)., (m2) > -I?^(„,t). 

Next we show how to normahze the problem at infinity. The use of the above sym- 
metries is necessary and it makes essential use of the second sheet of the Riemann 
surface (see also the Conclusion of this paper). 

Theorem 3.1. The function 

(3.21) m^{p)^ } w?{p,n,t) 

with m?{p,n,t) defined in (I3.13P is meromorphic away from S and satisfies: 
m\{p)^mUp)J^{p), peJ:, 

(3.22) (™?) > ("^2) > -%«,«)' 
m^{p*)^m^p)(^^ J) 

(3.23) TO^(oo+) = (1 *) , 
where the jump is given by 



r3 24) I^(vnt)-f l-^l' -Rip)e{p,n,t)e-'*(PA 

Setting R{z) = we clearly recover the purely periodic solution, as we should. 
Moreover, note 

(3.25) m^p)=(^±^ ^) + {^tt^ -i?+(n-l,0)^ + O(^). 

for p — (z, — ) near oo^. 

While existence of a solution follows by construction, uniqueness follows from 
Theorem El] and Remark EH 



Theorem 3.2. The solution of the Riemann-Hilbert problem of Theorem is 
unique. 



4. The stationary phase points and corresponding contour 

deformations 

The phase in the factorization problem p.lSp is t </> where cj) was defined in (|3.17p . 
Invoking (|2.15p and (|2.16p . we see that the stationary phase points are given by 

(4.1) ;Q(,„A,) + ^n(z-A,) = 0. 

j=o j=l 

Due to the normalization of our Abelian differentials, the numbers Aj, 1 < j < g, 
are real and different with precisely one lying in each spectral gap, say Xj in the 
j'th gap. Similarly, Aj, < j < g, are real and different and \j, 1 < j < g, sits in 
the j'th gap. However Ao can be anywhere (see [40l Sect. 13.5]). 

As a first step let us clarify the dependence of the stationary phase points on j . 
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Lemma 4.1. Denote by Zj{r/), < j < g, the stationary phase points, where rj = j. 
Set Ao = — oo and Xg+i — oo, then 

(4.2) \j < Zj{ri) < Aj+i 

and there is always at least one stationary phase point in the j 'th spectral gap. 
Moreover, Zj{r]) is monotone decreasing with 

(4.3) lim Zj{ri) = Xj+i and lim 2,(77) = ^j- 

Proof. Due to the normalization of the Abehan differential ^o+V^oo^ oo_ there is at 
least one stationary phase point in each gap and they are all different. Furthermore, 

^/ ^ lizj) ^ nLi(^j - ^fc) 

^ q'izj) + m'{zj) riLo.Mj '^i ~ ' 

where 

9 9 

q{z) = \{{z-\k), q{z) = \{{z-\k). 
k=0 k=l 

Since the points Afc are fixed points of this ordinary first order differential equation 
(note that the denominator cannot vanish since the Zj's are always different), the 
numbers Zj cannot cross these points. Combining the behavior as — > ±c» with 
the fact that there must always be at least one of them in each gap, we conclude 
that Zj must stay between Aj and Aj+i. This also shows z'^ < and thus Zj(ri) is 
monotone decreasing. □ 

In summary, the lemma tells us that we have the following picture: As j runs 
from —00 to +00 we start with Zg{ri) moving from 00 towards £'29+1 while the others 
stay in their spectral gaps until Zg{r]) has passed the first spectral band. After this 
has happened, Zg-i{r]) can leave its gap, while Zg{r]) remains there, traverses the 
next spectral band and so on. Until finally zo{r]) traverses the last spectral band 
and escapes to —00. 

So, depending on n/t there is at most one single stationary phase point belonging 
to the union of the bands a{Hq), say Zj{n/t). On the Riemann surface, there are 
two such points Zj and its flipping image z* which may (depending on n/t) lie in 
S. 

There are three possible cases. 

(i) One stationary phase point, say Zj, belongs to the interior of a band 
[E2j , E2j+i] and all other stationary phase points lie in open gaps. 

(ii) Zj — Zj — Ej for some j and all other stationary phase points lie in open 
gaps. 

(iii) No stationary phase point belongs to a{Hq). 
Case (i). Note that in this case 

(4.4) r{z,)/i= ^'='-\%^"'~"'K o. 

1-^23+2 (2j) 

Let us introduce the following "lens" contour near the band [i?2j, i?2j"+i] as shown 
in Figure [2j The oriented paths Cj = Qi U Cj2, C* = C*-^ U C*2 are meant to be 
close to the band [£'2^, i?2j+i]- 
We have 

Re((/)) > 0, in Dji, Re(9i) < 0, in Dj2. 
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+ 




+ 



Figure 2. The lens contour near a band containing a stationary 
phase point Zj and its flipping image containing z*. Views from 
the top and bottom sheet. Dotted curves lie in the bottom sheet. 



+ 




E2k E2k+1 



Figure 3. The lens contour near a band not including any sta- 
tionary phase point. Views from the top and bottom sheet. 

Indeed 

(4.5) lm{(j)') < 0, in [E2j,Zj], lm{(j)') > 0, in [zj,E2j+i] 

noting that is imaginary in [E2j-, E2j+i] and writing cp' = d(j)/dz. Using the 
Cauchy-Riemann equations we find that the above inequalities are true, as long as 
Cji,Cj2 are close enough to the band [E2j, E2j+i]. A similar picture appears in 
the lower sheet. 

Concerning the other bands, one simply constructs a "lens" contour near each of 
the other bands {E2k, E2k+i] and [-£^2fc' -^2fc+i] shown in Figure [S] The oriented 
paths Ck,C^ are meant to be close to the band [E2k, E2k+i]- The appropriate 
transformation is now obvious. Arguing as before, for all bands [E2k, E2k+i] we 
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will have 

Re(</.) < (>)0, in Dk, k > (<)j. 

Now observe that our jump condition (j3.24p has the following important factor- 
ization 

(4.6) J3 = {b-y^b+, 

where 

1 mc-*'l'\ h - ( ^ ^ 

1 / ' XRQe^'f' 1 

This is the right factorization for z > Zj{n/t). Similarly, we have 

where 

This is the right factorization for z < Zj{n/t). To get rid of the diagonal part 
we need to solve the corresponding scalar Riemann-Hilbert problem. Again we 
have to search for a meromorphic solution. This means that the poles of the scalar 
Riemann-Hilbert problem will be added to the resulting Riemann-Hilbert problem. 
On the other hand, a pole structure similar to the one of is crucial for uniqueness. 
We will address this problem by choosing the poles of the scalar problem in such a 
way that its zeros cancel the poles of m? . The right choice will turn out to be 
(that is, the Dirichlet divisor corresponding to the limiting lattice defined in (ILGI) ). 

Lemma 4.2. Define a divisor 2?p(„^t) of degree g via 

(4-8) ap^{'D0(n,t)) = ap„{'Dp.(n,t)) + S{n,t), 

where 

(4.9) Sein,t) = ^f log(l - |i?nc,. 

Jc{n/t) 

Then is nonspecial and ■K{vj{n,t)) — Vj{n,t) e M with precisely one in each 

spectral gap. 

Proof. Using (|2.15p one checks that 5i is real. Hence it follows from [JUl Lem. 9.1] 
that the Vj are real and that there is one in each gap. In particular, the divisor I?p 
is nonspecial by j40[ Lem. A. 20]. □ 

Now we can formulate the scalar Riemann-Hilbert problem required to eliminate 
the diagonal part in the factorization <\A.7\ : 

d+{p,n,t) = d^{p,n,t){l-\R{p)\^), peC{n/t), 

(4.10) (rf) >-%„,*), 
d(oo_|-, n, t) = 1, 

where C{n/t) = E n 7r~^((— oo, Zj{n/t)). Since the index of the (regularized) jump 
is zero (see remark below), there will be no solution in general unless we admit g 
additional poles (see e.g. [351 Thm. 5.2]). 
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Theorem 4.3. The unique solution of (|4.10|) is given by 




e{z{n,t) + S{n,t)) 9{z{p,n,t)) 



X 



9{z{n,t)) e{z{p,n,t)+S{n,t)) 



X exp — / log(l - \Rf)ujpoo 

WtTI Jc(n/t) 



) 



where S_{n, t) is defined in (j4.9p and ujpq is the Abelian differential of the third kind 
with poles at p and q. 

The function d(p) is meromorphic m M \ E with first order poles at i'j{n,t) and 
first order zeros at fij{n,t). Also d{p) is uniformly hounded in n,t away from the 
poles. 

In addition, we have d{p) — d{p). 

Note that this formula is different (in fact much simpler) from the explicit solu- 
tion formula from Rodin [35l Sec. 1.8]. It is the core of our explicit formula (|1.6|) 
for the limiting lattice. 

Proof. On the Riemann sphere, a scalar Riemann-Hilbert problem is solved by 
the Plemelj-Sokhotsky formula. On our Riemann surface we need to replace the 
Cauchy kernel by the Abelian differential of the third kind a;poo+- But now 
it is important to observe that this differential is not single-valued with respect 
to p. In fact, if we move p across the ai cycle, the normalization J^^ ujpac^ — 
enforces a jump by 2TriQ. One way of compensating for these jumps is by adding to 
Wpoo-i- suitable integrals of Abelian differentials of the second kind (cf. [351 Sec 1.4] 
respectively Section Since this will produce essential singularities after taking 
exponentials we prefer to rather leave u;poo+ as it is and compensate for the jumps 
(after taking exponentials) by proper use of Riemann theta functions. 
To this end recall that the Riemann theta function satisfies 



where r is the matrix of 6-periods defined in (|2.7|) and (., ..) denotes the scalar 
product in (cf., e.g. [15] or [40l App. A]). By definition both the theta functions 
(as functions on M) and the exponential term are only defined on the "funda- 
mental polygon" M of M and do not extend to single-valued functions on M in 
general. However, multi-valuedness apart, d is a (locally) holomorphic solution of 
our Riemann-Hilbert problem which is one at oo+ by our choice of the second 
pole of the Cauchy kernel LOp oo+. The ratio of theta functions is, again apart from 
multi-valuedness, meromorphic with simple zeros at fij and simple poles at i)j by 
Riemann's vanishing theorem. Moreover, the normalization is chosen again such 
that the ratio of theta functions is one at oo+. Hence it remains to verify that 
(|4.1ip gives rise to a single-valued function on M. 

Let us start by looking at the values from the left /right on the cycle b^. Since our 
path of integration in z_{p) is forced to stay in M, the difference between the limits 
from the right and left is the value of the integral along a^. So by (|4.12p the limits 
of the theta functions match. Similarly, since ujp oc^ is normalized along at cycles, 
the limits from the left/right of ujpoo^ coincide. So the limits of the exponential 
terms from different sides of be match as well. 



(4.12) 
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Next, let us compare the values from the left/right on the cycle ag. Since our 
path of integration in z_{p) is forced to stay in M, the difference between the limits 
from the right and left is the value of the integral along be. So by (|4.12|) the limits 
of the theta functions will differ by a multiplicative factor exp{2TTi6e). On the other 
hand, since Wpoo+ is normalized along ae cycles, the values from the right and left 
will differ by —2TTiQ. By our definition of S_ in (j4.9p . the jumps of the ration of 
theta functions and the exponential term compensate each other which shows that 
(|4.11|) is single-valued. 

To see uniqueness let d be a second solution and consider d/d. Then d/d has no 
jump and the Schwarz reflection principle implies that it extends to a meromorphic 
function on M. Since the poles of d cancel the poles of d, its divisor satisfies 
(d/d) > —Dp,. But 2?^ is nonspecial and thus d/d must be constant by the Riemann- 

Roch theorem. Setting p = oo^ we see that this constant is one, that is, d = d as 
claimed. 

Finally, d{p) — d{p) follows from uniqueness since both functions solve ()4.10|) . □ 

Remark 4.4. Once the last stationary phase point has left the spectrum, that is, 
once C{n/t) = Yj, we have d{'j 
A = A+{n,t)A^{n,t) = r(oo). 



once C{n/t) = Yi, we have d{p) = A ^T{z)^^, p — {z,zt) (compare [42] ). Here 



In particular. 



(4.13) 



J{z{n-l,t)) e{z{n,t)+6{n,t)) 
a^c50_,n,rj 0(^z{n,t)) e{z{n - l,t) + S{n,t)) 

X exp ( / log(l - |i?|^)woo_ c 

Jc(n/t) 



since _z(oo_, n, t) = z(oo+, ri— 1, t) — z(ri— 1, t). Note that (i(oo_, n, t) = d{oo- , n, t) - 
d(oo_,n,t) shows that d{oo-,n,t) is real- valued. Using (|2.15p one can even show 
that it is positive. 

The next lemma characterizes the singularities of d{p) near the stationary phase 
points and the band edges. 

Lemma 4.5. For p near a stationary phase point Zj or z* (not equal to a band 
edge) we have 

(4.14) d(p) = (z-z,)±'''e±(z), p={z,±), 
where e^(z) has continuous limits near zj and 

(4.15) ^=_i.iog(i_|i?(z^.)|2)>o. 

Here {z — zj)^^'^ = exp(±ij/log(z — zj)), where the branch cut of the logarithm is 
along the negative real axis. 

For p near a band edge G C'{n/t) we have 

(4.16) d(p) =r±i(z)e~±(z), p^{z,±), 

where e^(z) is holomorphic near Ek if none of the Vj is equal to Ek and e±{z) has 
a first order pole at Ek = Vj else. 
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Proof. The first claim we first rewrite (|4.11|) as 

\ 0{z{n,t)+S{n,t)) 0{z{p,n,t)) 



d{p, n, t) — exp iv I 



Next observe 

(4.18) If ujpp, ^±log{z- zj)±a{zj) + 0{z- zj), p^{z,±), 

^ ■IC{n/t) 



where a{zj) G M, and hence 
(4.19) 

/ UJp 00+ ^ ±^Og{z- Zj)±a{Zj) + - UJoo^oo++0{z-Zj), p={z,±), 

JC(n/t) ^ Jc(n/t) 

from which the first claim follows. 
For the second claim note that 

r(oo) \r(z)"i, p = (z,-)en_, 

satisfies the (holomorphic) Riemann-Hilbert problem 

i+(p) =t_(p)(i-|i?(p)n, pes, 

t[QO + ) = 1. 

Hence d{p)/t{p) has no jump along C{n,t) and is thus holomorphic near C{n/t) 
away from band edges Ek = Vj (where there is a simple pole) by the Schwarz 
reflection principle. □ 

Furthermore, 
Lemma 4.6. We have 



(4.20) e±(z) = eT(z), p = (z, ±) e I]\C(n/t), 

and 



e^(2^i) =exp {wa{zj) + y ^ 



C(n/t) ) 

0{z{n, t) + 5_{n, t)) 0{z{zj,n, t)) 
^ 0iz{n,t)) d{z{zj,n,t) + S{n,t)) ^ 



+ Wo 



where 

(4.22) a(zj)^ lini i / Wpp. - log(7r(p) - z^). 

2 Jc(n/t) 

77ere a(zj) G R anrf cjpp. is rea/ whereas UJoc- oo+ is purely imaginary on C{n/t). 
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Proof. The first claim follows since d{p*) — d{p) — d{p) for p e T,\C{n/t). The 
second claim follows from (|4.17p using Jc{n/t) f ^p°°+ ^ 5 Ic{n/t) f (^pp+^oo. oo+) 
for symmetric fmictions f{q) = f{q*)- O 

Having solved the scalar problem above for d we can introduce the new Riemann- 
Hilbert problem 

(4.23) m*{p) = d{oo^)-'m'ip)Dip), Dip) ^ (^^^P*^ f^^^ 

where d*{p) = d{p*) is the unique solution of 

dXip) - d*_{p){l ~ \R{p)\')-\ P e C{n/t), 

d*(oo_) = 1. 

Note that 



det(D(p)) = dip)d{p*) = d((x)_) n 



g 



Z — V 



Then a straightforward calculation shows that rr& satisfies 

= fnt{p)J'^{p), pes, 

(4.24) (mt) > -%n,t)-, ("^2) 
m*(p*) = m^(p) 

m*(oo+) (1 *) , 
where the jump is given by 

(4.25) J*ip) = D_{p)-'j^p)D+{p), p€E. 

In particular, has its poles shifted from jlj{n,t) to Vj{n,t). 
Furthermore, can be factorized as 

(4.26) J'^{l~JS<b ~^^^''~'')={b-)-'h, peJ:\C(n/t), 



^^RQe*'!' 1 
where b± — D^^b±D, that is, 

'1 SrR&e^'' 



(4.27) b- 
for tt{p) > Zj{n/t) and 



1 }' V^i?ee*'> 1 



(4-28) ^'=L*_.. - U(i?-)-^S+, peC(n/t), 



where i3± = Dj.^B±D±, that is. 



/ 1 0\ - /i -it-Re e-*-^^ 



for 7r(p) < Zj{n/t) 
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Note that by d{p) — d{p) we have 



respectively 



oiN d+{p) d+{p) 1 dl.ip) 

^'■''^ d^p) - mr^mw " ' ""^"Z'^- 

We finally define by 

= m^B^'^, pe Dk, k < j, 

rrt' = m'^BZ^, p G Dl, k < j, 

m^ = m'^B^^, peDji, 

m^^m^BZ\ peD*^, 

(4.32) m'^ = m%\ p e Dj2, 

m'^ = m'^bz\ peD*2, 

— rn^b^^ , p G Dk, k > j, 

= m'^bZ^, p e Dl, k > j, 

rrt' = m^, otherwise, 

where we assume that the deformed contour is sufficiently close to the original one. 
The new jump matrix is given by 

n, t) — m^_{p, n, t)J^{p, n, t), 

J^^B+, peCk, k<j, 

J^-SZ\ peCt, k<j, 

= B+, peCji, 

(4.33) J^ = SZ\ peC*!, 

J^ = 6+, peCj2, 
.f = bz\ peC*2, 

^b+, pe Ck, k > j, 

j^^bz\ peCl k>j. 

Here we have assumed that the function R{p) admits an analytic extension in the 
corresponding regions. Of course this is not true in general, but we can always 
evade this obstacle by approximating R{p) by analytic functions in the spirit of |6j. 
We will provide the details in Section [6l 

The crucial observation now is that the jumps on the oriented paths Ck, 
are of the form I + exponentially small asymptotically as i — > oo, at least away 
from the stationary phase points zj , z* . We thus hope we can simply replace these 
jumps by the identity matrix (asymptotically as t ^ oo) implying that the solution 
should asymptotically be given by the constant vector (l l) . That this can in fact 
be done will be shown in the next section by explicitly computing the contribution 
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of the stationary phase points thereby showing that they are of the order 0{t ^/^), 
that is, 



(p) = (i i) + orV2) 



uniformly for p a way from the jump contour. Hence ah which remains to be done 
to prove Theorem 11.11 and Theorem 11.31 is to trace back the definitions of m* and 
and comparing with p.25p . First of ah, since and coincide near cx)_ we 
have 

m\p)^{l l)+0(<-i/2) 

uniformly for p in a neighborhood of cxd_. Consequently, by the definition of to^ 
from ()4.23p . we have 

m^p) = d(oo_) {d{p*)'^ d{py^) + 0{t-'/^) 

again uniformly for p in a neighborhood of cx)_ . Finally, comparing this last identity 
with p.25|) shows 

(4.34) A+{n,t)^ = d{oo-,n,t) + 0(t-^/'^), B+{n,t) = -di{n,t) + 0{t-^/^), 
where di is defined via 

d{p) = 1 + — + 0{\), p = {z, +) near cx3+ 
z z^ 

Hence it remains to compute c?i. Proceeding as in [40, Thm. 9.4] respectively [HI 
Sec. 4] one obtains 



di=-^ [ iog(i - 

^TTl Jc{n/t) 

1 d ^^^fO{z{n,s)+S{n,t)) 



2ds V S{zin,s)) 
where fio is the Abelian differential of the second kind defined in (|2.16p . 



Case (ii). In the special case where the two stationary phase points coincide (so 
Zj = z* = Ek for some k) the Riemann-Hilbert problem arising above is of a 
different nature, even in the simpler non-generic case \R{Ek)\ < 1. In analogy to 
the case of the free lattice one might expect different local asymptotics expressed 
in terms of (generalized) Painleve functions. In the case \R{Ek)\ < 1 the two 
crosses coalesce and the discussion of Section |B] goes through virtually unaltered. If 
\R{Ek ) I = 1 the problem is singular in an essential way and we expect a " collisionless 
shock" phenomenon in the region where Zj{n/t) ~ E^, similar to the one studied 
in [T], [3], [53]. The main difficulty arises from the singularity of jiq^- We expect 
that an appropriate "local" Riemann-Hilbert problem however is still explicitly 
solvable. It is definitely interesting to understand the actual contribution of the 
band edges and the way it is different from the free case due to the fact that the 
underlying Riemann surface has nonzero genus. We plan to investigate this issue in 
a future publication. But in the present work, we will assume that the stationary 
phase points stay away from the Ek ■ 
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Case (iii). In the case where no stationary phase points he in the spectrum the 
situation is similar to the case (i). In fact, it is much simpler since there is no 
contribution from the stationary phase points: There is a gap (the j-th gap, say) 
in which two stationary phase points exist. We construct "lens- type" contours Ck 
around every single band lying to the left of the j-th gap and make use of the 
factorization — {b^)~^b+. We also construct "lens- type" contours Ck around 
every single band lying to the right of the j-th gap and make use of the factorization 
J3 = {B^y^B+. Indeed, in place of (|432t we set 

= m'^B^^, p e Dk, k < j, 

= m'^BZ^, p e D'l, k < j, 

(4.35) TO5 = m^6;\ peDk,k>j, 

rrt'^m'^hZ^, p e D*^, k > j, 



'k^ 

m" — rrC 



n5 _ ^4 otherwise. 



It is now easy to check that in both cases (i) and (iii) formula (j4.13p is still true. 

5. The "local" Riemann-Hilbert problems on the small crosses 

In the previous section we have shown how the long-time asymptotics can be 
read off from the Riemann-Hilbert problem 

?Ti+(j', n, t) — mt{p, n, t)J^{p, n, t), p £ E^, 
{ml) > -V^(^n,t)' , {m2) > -'Do{n,t), 
{p* ,n,t) = nt' {p,n,t) 
(5.1) m5(oo+,n,i) (1 *) . 

In this section we are interested in the actual asymptotic rate at which m^{p) — )■ 
(l l). We have already seen in the previous section that the jumps on the 
oriented paths Ck, C^ for k ^ j are of the form I -I- exponentially small asymptot- 
ically as i — >■ oo. The same is true for the oriented paths Cji,Cj2,C*i,C*2 at least 
away from the stationary phase points Zj, z*. On these paths, and in particular 
near the stationary phase points (see Figure 14]), the jumps read 



1 - 



d R'e' 



= B_^ ^ ( t ^ j, peC*j, 

1 

i^ee*"^ 1 



J = 6+ = ( ^. , peCj2, 



(5.2) ^' = ^=' = (o 1 p^c;2. 

Note that near the stationary phase points the jumps are given by (cf. Lemma 

X 2ii/ 




1 
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Figure 4. The small cross containing the stationary phase point 
Zj and its flipping image containing z*. Views from the top and 
bottom sheet. Dotted curves lie in the bottom sheet. 




where (cf. ((HT^ and (liT^ ) 



(5.4) r = i?(z^.)e(z„n,t)^ ^IMJ . 

The error terms will satisfy appropriate Holder estimates, that is 

(5.5) \\B+{p)~B+{p)\\<C\z-z,r, 

for any a < 1 and similarly for the other matrices. 

To reduce our Ricmann-Hilbert problem to the one corresponding to the two 
crosses we proceed as follows: We take a small disc D around Zj{n/t) and project it 
to the complex plane using the canonical projection tt. Now consider the (holomor- 
phic) Riemann-Hilbert problem in the complex plane with the very jump obtained 
by projection and normalize it to be I near oo. Then, as is shown in [6] (see |281 
Thm. A.l]), the solution is of the form 

(5.6) M(z)=I+^-L + o(i-), 
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where 



\m 

(5.7) f3{t) — yfl}e^'''^^'^~^'^^^^^^^'^^^^^'^'^^^^e~^*'^^^^h~^'^ . 

Now we lift this solution back to the small disc on our Riemann-surface by setting 
M (p) = M{z) ioi pe D and M{p) = M{z) for pG D*. Then 

(5.8) n.^ip)^hy''''^'^' 

Ito°(p), else. 

has no jump inside D U D* but jumps on the boundary given by 

(5.9) mlip) :^mtip)Ar\p), pedDUdD*. 

The remaining jumps are unchanged. In summary, all jumps outside D U dD* . are 
of the form I + exponentially small and the jump on dD U dD* is of the form 

I+0(t-l/2). 

In order to identify the leading behaviour it remains to rewrite the Riemann- 
Hilbert problem for as a singular integral equation following Section [X] Let the 
operator C^e : L'^{T,^) -> L2(S^) be defined by 

(5.10) a,./ = C_(/z(;6) 
for a vector valued /, where — — I and 

If- ■ fn-'°°+ n \ 



(5.11) (C±/)(g)= lim^,:^^/ /ilf, = 



r!f'°°- y ' 



are the Cauchy operators for our Riemann surface. In particular, Vlj; '^ is the Cauchy 
kernel given by 

(5.12) f]f'«=c.p, + ^/f''(p)C„ 

i=i 

where 

rp 



(5.13) if\p)=Y.^AiL) r 



Here w^^o is the (normalized) Abelian differential of the second kind with a second 
order pole at q (cf. Remark [Ol below) . Note that lf''^{p) has first order poles at 
the points £;_. 

The constants Cji{v) are chosen such that fi^''' is single valued, that is, 

-1 



(5.14) (c»(^:))i<^_fe<g = ^ck{i)-T 



J-1 



where Ck{j) are defined in (|2.6p (cf. Lemma [A. 3[) . 

Remark 5.1. The Abelian differential Upq is explicitly given by 

_ ( <12+<12(P) _ <+2+<+2(g) , \ _±_ 

""'''y 2(.-.(p)) 2in-n{q)) + ^'^"V<+2' 



24 



S. KAMVISSIS AND G. TESCHL 



where Ppq{z) is a polynomial of degree g — I which has to be determined from the 
normalization J^^ ujpp* =0. For q = oo± we have 



V 2(7r - Trip)) 2 rIH. 



Consider the solution fj,^ of the singular integral equation 

(5.15) M=(l l)+C^6M in ^^(E*^). 
Then the solution of our Riemann-Hilbert problem is given by 

(5.16) m6(p) = (l 1) + ^ 

By Ijw^lloo = 0{t~^^^) Neumann's formula implies 

(5.17) M'(3) = (I-C^a)-i(l 1) = (1 l)+0(f-V2). 
Moreover, 

f- — 7172 +0(t-"), pedD, 

(5.18) w^(p) = l y h P 



S6 



Z-Zi tl/2 



Hence we obtain 



m 



(I l)M. 1 r 

JdD' - Zj 



.1 1 
fl 1 



tV2 27ri 

) ^175 — ^175 — ) + ) 



(5.19) 



+ 0(i-"). 



^ (i/3nr'°°"(^i)-i/3i^f'°°+(^l) -i^J^f'°°-(^,) + i/34'°°-(^;)) 



Note that the right-hand side is real- valued for p e tt ^(R)\S since Q!^'°°^{q) 
n='°°='=(g) implies 



(5.20) f2i'°°±(4) =0f'°°*(^,), pe7r-i(M)\E. 
Since we need the asymptotic expansions around oo_ we note 
Lemma 5.2. We have 

(5.21) ^,oo^^^.)=^^+A^l + o{^) 

for p = (z, — ) near cx)_j where 

/oo_ 
c^p„oa(^,) 
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Af = OJoo-,o{^j) + ^Cfc£(i>)'^!/j,o(oo_)Cfe(2j) 
kj 

(5.23) = Woo_,o(zj) -^Cke{E)^i>*.,o{^+)Ck{zj). 

k,e 

Proof. To see n^_{zj) = Q.^_{zj) note Ckiiv*) = -cu{m) and /^~W;y;,o = 

Observe that since CkpXM) G and uJa^^o G K we have G iM. 

Remark 5.3. TVote that the Abelian integral appearing in the previous lemma is 
explicitly given by 

-7r^+^ + \ E-iV E,n^ + P^-A-) + <% , 

(5.24) t^oo_,o = Y/2 "-^^ 

with Poo_.o polynomial of degree g — I which has to be determined from the nor- 
malization. 
Similarly, 

<+2 + <+2('>) + - ^) + Pc'A^) ■ - 

(5.25) ^,,0 ^ "-^^'^ dn, 

2(^-z.)2i?i/^2 

with Pofl a polynomial of degree g—1 which has to be determined from the normal- 
ization. 

As in the previous section, the asymptotics can be read off by using 

(5.26) m3(p) = d(oo_)m6(p) [ S(ry M 

V " dip)] 

for p near oo- and comparing with p.25p . We obtain 



(5.27) ^An,tf-^A^ + J^, (iM - i/3Af ) + Oit^) 



and 



(5.28) B+in,t) = -di - (i/3Af -i/3Af ) +0(i-"). 

6. Analytic Approximation 

In this section we want to show how to get rid of the analyticity assumption 
on the reflection coefficient R{p). To this end we wih spht R{p) into an analytic 
part Ra^t plus a small rest i?^,* following the ideas of |B| (see also [28j Sect. 6]). 
The analytic part will be moved to regions of the Riemann surface while the rest 
remains on E = 7r~^ (cr(iJq)) . This needs to be done in such a way that the rest is 
of 0{t~^) and the growth of the analytic part can be controlled by the decay of the 
phase. 
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In order to avoid problems when one of the poles Uj hits S, we have to make the 
approximation in such a way that the nonanalytic rest vanishes at the band edges. 
That is, split R according to 

R^)^R{E,,) '-\ +R{E2,+,) 



i?2j + l — E2j E2j — E2j + l 



(6.1) ±^z-E2j^/z-E2j+iRip), P={z,±), 

and approximate R. Note that if i? e C'(E), then R e C'^^'^iY.). 

We will use different splittings for different bands depending on whether the 
band contains our stationary phase point Zj[n/t) or not. We will begin with some 
preparatory lemmas. 

For the bands containing no stationary phase points we will use a splitting based 
on the following Fourier transform associated with the background operator Hq. 
Given i? G C'(S) we can write 

(6.2) i?(p) = ^i?(ri)V',(p,n,0), 

where V'g(Pj x, t) denotes the time-dependent Baker-Akhiezer function and (cf. [lOj . 

ini) 



1 f in?-i(^(p)-Aij) 

If we make use of p.l2p . the above expression for R{p) is of the form 

(6.4) R{p) ^^R{n)eq{p,n,0)exp{ink{p)). 

where k{p) = — i /^^ Wcx)+ oo_ and 9q{p,n,t) collects the remaining parts in ()2.12p . 
Using k{p) as a new coordinate and performing / integration by parts one obtains 

, A/ s , const 

(6.5) |i?(n)| < 



1 + H 
provided i?e C'(S). 

Lemma 6.1. Suppose R £ ^^(Z), n'-R{n) e f{Z) and let (3 > be given. Then we 
can split R{p) according to 

R{p)=RaM+RrAp)^ 

such that Ra^tip) 'is analytic for in the region < lm{k{p)) < e and 

(6.6) \RaAp)e'^'\ = Oit-'), < Im(fc(p)) < e, 

(6.7) \RrAp)\ = 0{t-'), pes. 

Proof. We choose 

OO 

Ra.tip) = ^ R{n)0q{p,n,O)exp{mk{p)) 
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with N(t) = l^t\ for some positive /3o < /3- Then, for < lm{k{p)) < e, 

oo 

|-Ra,t(fc)e-''*| < Ce-^* J2 |i?(ri)|e-i'"('=(''»" 

n=-N{t) 

<Ce-'^*e^(*)-n|F||i = ||i?||ie-(^-'3«)*, 
which proves the first claim. Similarly, for p G S, 



n=JV(t) + l 



□ 



For the band which contains Zj{n/t) we need to take the small vicinities of the 
stationary phase points into account. Since the phase is cubic near these points, we 
cannot use it to dominate the exponential growth of the analytic part away from E. 
Hence we will take the phase as a new variable and use the Fourier transform with 
respect to this new variable. Since this change of coordinates is singular near the 
stationary phase points, there is a price we have to pay, namely, requiring additional 
smoothness for R{p). 

Without loss of generality we will choose the path of integration in our phase 
4>{p)^ defined in p.l7p . such that (l){p) is continuous (and thus analytic) in Dj^i 
with continuous limits on the boundary (cf. Figure [2]). We begin with 

Lemma 6.2. Suppose R{p) G C'^(S). Then we can split R{p) according to 

(6.8) R{p) = Ro{p) + {n{p)-n{zj))H{p), peSn^i, 

where Ro{p) is a real rational function on M such that H{p) vanishes at Zj, z* of 
order three and has a Fourier series 

97ri 

(6.9) Hip) . |: ^(n)e-.(.), = ^(^^^^ > 0^ 
with nH{n) summahle. Here (j) denotes the phase defined in ()3.17p . 

1 /2 

Proof. We begin by choosing a rational function Ro{p) — a{z) + 6(2:)i?2g+2(?') '^ith 
p = {z,±) such that a(z), b{z) are real- valued polynomials which are chosen such 
that a{z) matches the values of Re(i?(p)) and its first four derivatives at Zj and 
i^^6(z)i?2g+2(?') matches the values of Im(i?(p)) and its first four derivatives at Zj. 
Since R{p) is we infer that H{p) S C'*(S) and it vanishes together with its first 
three derivatives at Zj, z*. 

Note that (j){p)/i, where (j) is defined in p.l7p has a maximum at z* and a 
minimum at Zj. Thus the phase (j){p)/i restricted to S H Dj i gives a one to one 
coordinate transform S H Dj^i — > [0(z*)/i, 0(zj)/i] and we can hence express H{p) 
in this new coordinate. The coordinate transform locally looks like a cube root near 
Zj and z* , however, due to our assumption that H vanishes there, H is still in 
this new coordinate and the Fourier transform with respect to this new coordinates 
exists and has the required properties. □ 



Moreover, as in Lemma l6. II we obtain: 
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Lemma 6.3. Let H{p) be as in the previous lemma. Then we can split H{p) 
according to H{p) — Ha.t{p) + Hr.t{p) such that Ha.t{p) is analytic in the region 
Re((/)(p)) < and 

(6.10) |i?a,t(p)e^(p)*/2|=0(l), pe:D~T, \HrM\^0{t-^),pe^. 

Proof. We choose Ha,t{p) = Yln=-K(t) H {n)e"'^o'f'^P'> with K{t) = lt/{2ujo)\. Then 
we can proceed as in Lemma l6. II 

|i?a,t(p)e'^(P^*/^| < ||F||i|e--^(*)'^''"^(P)+'^(f)*/2| < ||ff||i 

and 



\Hr,tip)\ < 



1 



Kit) 



J2 n\H{-n)\ < 



C 



n=K{t) + l 



□ 



Clearly an analogous splitting exists for p € S n Dj2 . 

Now we are ready for our analytic approximation step. First of all recaU that our 
jump is given in terms 6± and B± defined in (I4.27P and (|4.29p . respectively. While 
h± are already in the correct form for our purpose, this is not true for B± since they 
contain the non-analytic expression |r(p)p. To remedy this we will rewrite B± in 
terms of the left rather than the right scattering data. For this purpose let us use 
the notation Rr{p) = R+{p) for the right and Ri{p) = R-{p) for the left reflection 
coefficient. Moreover, let dr{p,x,t) = d{p,x,t) and di{p,x,t) = T{p)/d{p,x,t). 

With this notation we have 



(6.11) 
where 



Ap) 



h-{p)-^b+{p), 
B-{p)-^B+{p), 



Ti{p) > Zjin/t), 
tt{p) < Zj{n/t), 



6_ = 



^^R^ip*)eip*)e-*^(r>^ 



and 



i^^i?.(p)e(rte-*^W 



B_ = 

Using p.7p we can write 

B^ = 

B, = 



-(p',xS) Rr(p)<Si{p) „tMv) 1 ' 

dr,+ (p,x,t) R^ip*)e{p') td>(v)'' 
d^,+ {p',x,t) |T(p)P ^ 



1 il^MP*MP*)e''^''^' 

1 



Now we split Rr{p) — Ra,t{p)+Rr,t{p) by splitting Rr{p) defined via (|6.ip according 
to Lemma l6.ll for 7r(p) e [i?2fc, ^'2fc+i] with k < j (i.e., not containing Zj{n/t)) 



ASYMPTOTIC STABILITY OF THE PERIODIC TODA LATTICE 



29 



and according to Lemma 16.31 for 7r(p) G [E2j , Zj{n/t)]. In the same way we split 
Ri{p) = Ra,t{p) + Rr,t{p) for ^b) £ [^j('^A))-S2j+i] and tt{p) e [i?2fc, -^2/0+1] with 
k > j. For /3 in Lemma 16.11 we can choose 



(6.12) 



minpeCfc -R-e((?!)(p)) > 0, 
minpgCfc Re(0(p)) > 0, 



tt{p) > Zj{n/t), 
tt{p) < Zj{n/t). 



In this way we obtain 

B±ip) = Ba,u±iP)Br,t,±iP) 



br,t,±{p)ha,t,±{p), 
= Br,t,±{p)Ba,t,±{p)- 



Here &a,t,±(p), W,t.±{p) (resp. Ba,t,±(jp), Br^t,±{p)) denote the matrices obtained 
from b±{p) (resp. S±(p)) by replacing Rrip) (resp. i?i(p)) with Ra,t{p), Rr^ip), 
respectively. Now we can move the analytic parts into regions of the Riemann 
surface as in Section |4] while leaving the rest on E. Hence, rather than (|4.33p . the 
jump now reads 



(6.13) 



J'ip) = 



ba,t,+ {p), 

Kt,-{pr\ 
br^t-{py%^t.+{p), 

Ba,t.+ {P), 

Ba,tApr\ 

Br^t,-{Py^Br^t,+ 



ip), 



P e Cfe, 

P&C*„ 

pes, 

p e Cfe, 

pec,*, 

pes. 



7r(p) > Zj{n/t), 
tt{p) > Zj{n/t), 

7r{p) > Zj{n/t), 
tt{p) < Zj{n/t), 
Tr{p) < Zj{n/t), 

Tr{p) < Zj{n/t). 



By construction Ra.t{p) = Ra{p) + i^^ip) ~ '^{zj))Ha.t{p) will satisfy the required 
Lipschitz estimate in a vicinity of the stationary phase points (uniformly in t) and 
the jump will be J^{p) = I + 0{t^^). The remaining parts of E can be handled 
analogously and hence we can proceed as in Section [5l 



7. Conclusion 

We have considered here the stability problem for the periodic Toda lattice under 
a short-range perturbation. We have discovered that a nonlinear stationary phase 
method (cf. [6], [22]) is applicable and as a result we have shown that the long- 
time behavior of the perturbed lattice is described by a modulated lattice which 
undergoes a continuous phase transition (in the Jacobian variety). 

We have extended the well-known nonlinear stationary phase method of Deift 
and Zhou to Riemann-Hilbert problems living on a two-sheeted Riemann surface. 
We were forced to tackle such Riemann-Hilbert problems by the very problem, since 
there is no way we could use the symmetries needed to normalize the Riemann- 
Hilbert problem of Section [3] without including a second sheet. We believe that 
the fact that our analysis involves Riemann-Hilbert problems living in a Riemann 
surface is the main novelty of our contribution. 

Although the most celebrated applications of the deformation method initiated 
by [5] for the asymptotic evaluation of solutions of Riemann-Hilbert factorization 
problems have been in statistical mechanics (via orthogonal polynomials, random 
matrices and combinatorial probability), most mathematical innovations have ap- 
peared in the study of nonlinear dispersive PDEs or systems of ODEs (cf. [6], [9], 
p5] ) . It is thus interesting that another mathematical extension of the theory (set- 
ting the factorization problem on a Riemann surface rather than the complex plane) 
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arises in the study of an innocent looking stability problem for the periodic Toda 
lattice. 

Let us also remark here that there are other problems in the literature involving 
Riemann-Hilbert problems living in a Riemann surface (see for example [4]). We 
thus expect our methods to have a wide applicability. 



In the complex plane, the solution of a Riemann-Hilbert problem can be reduced 
to the solution of a singular integral equation (see [2]). In our case the underlying 
space is a Riemann surface M. The purpose of this appendix is to generalize this 
approach to Riemann-Hilbert problems of the type 



Concerning the jump contour E and the jump matrix J we will make the following 
assumptions: 

Hypothesis H. A.l. Let S consist of a finite number of smooth oriented finite 
curves in M which intersect at most finitely many times with all intersections being 
transversal. The divisor 2?^ is nonspecial. The contour E does neither contain 
oo± nor any of the points jl and that the jump matrix J is nonsingular and can be 
factorized according to J = bZ^b^ = (I — ii;_)^"'^(I + w+), where w± — i{b± — I) 
are continuous. 

Remark A. 2. (i). We dropped our symmetry requirement 



here since it only is important in the presence of solitons. However, if both S and 
w± are compatible with this symmetry, then one can restrict all operators below to 
the corresponding symmetric subspaces implying a symmetric solution. Details will 
be given in '28' . 

(a). The assumption that that none of the poles fi lie on our contour E can 
be made without loss of generality if the jump is analytic since we can move the 
contour a little without changing the value at oo^ (which is the only value we are 
eventually interested in). Alternatively, the case where one (or more) of the poles 
pLj lies on E can be included if one assumes that w± has a first order zero at fij . 
In fact, in this case one can replace fi{s) by /i(s) = (7r(s) — fij)fi(s) and w±{s) by 
w±{s) = (7r(s) - Hjy^w±{s). 

Otherwise one could also assume that the matrices w± are Holder continuous and 
vanish at such points. Then one can work with the weighted measure — ii?2g+2(?')'^^ 
on T,. In fact, one can show that the Cauchy operators are still bounded in this 
weighted Hilbert space (cf. |17l Thm. 4.1]j. 

Our first step is to replace the classical Cauchy kernel by a "generalized" Cauchy 
kernel appropriate to our Riemann surface. In order to get a single valued kernel 
we need again to admit g poles. We follow the construction from [351 Sec. 4]. 



Appendix A. A singular integral equation 



(A.1) 



m+{p) m^{p)J{p), _p e E, 
(mi)>-X'p,., {m2)>-T>p_, 
m(oo+) ~ mo £ C^. 



(A.2) 
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Lemma A. 3. Let Vp, be nonspecial and introduce the differential 



3 

(A.3) nf- = cjpoo+ + ^ lf{p)Cj, 

where 

3 rp 

Here UJq^Q is the (normalized) Abelian differential of the second kind with a second 
order pole at q (cf. Remark \5.3\) and the matrix Cji is defined as the inverse matrix 
of r]f,{jlj), where Q = r]i(z)dz is the chart expression in a local chart near flj (the 
same chart used to define Uj^.^a)- 

Then is single valued as a function of p with first order poles at the points (i. 

Proof. Note that I^{p) has first order poles at the points ji hence it remains to 

show that the constants Cjt{fl) are chosen such that is single valued (cf. the 
discussion in the proof of Theorem 14.31) . That is, 

/3 r 3 

dif = ^ Cji / ujf^^fi = ^ CjirikiM = Sjk, 

where = Vk{z)dz is the chart expression in a local chart near fii (here the bk 
periods are evaluated using the usual bilinear relations, see [T5j Sect. III. 3] or [40l 
Sect. A. 2]). That the matrix r]k{fii) is indeed invertible can be seen as follows: 
If X]fc=i Vkifi-p^Ck = for 1 < ^ < g, then the divisor of C = X]fc=i ^fcCfe satisfies 
(C) ^ T^fi- But since we assumed the divisor Vf^ to be nonspecial, i{T>fi) — 0, we 
have C = implying Cfc = 0. □ 

Next we show that the Cauchy kernel introduced in (|A.3|) has indeed the correct 
properties. We will abbreviate LP{Y.) = LP{T,,C'^)- 

Theorem A.4. Set 

(A.5) 

and define the matrix operators as follows. Given a 2 x 2 matrix f defined on S 
with Holder continuous entries, let 

(A.6) (Cf)ip)^^ [ fni for p^E, 

and 

(A.7) (C±/)(g)= lim (C/)(p) 

from the left and right ofT, respectively (with respect to its orientation). Then 
(i) The operators C'± are given by the Plemelj formulas 
{C+f){q)-iC^f){q)^f{q), 

{C+f){q) + {C^f){q) = --f fni 

TTl . -p 
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and extend to bounded operators on Here f denotes the principal 

value integral, as usual. 

(ii) Cf is a meromorphic function off S, with divisor given by {{Cf)ji) > 
-P^. and {{C,f)j2) > ~V^. 

(iii) (C/)(^+)-0. 

Proof. In a chart z = z{p) near go G S, the differential Vlj ~ {j^^q) + 0{l))dz 
and hence the first part follows as in the Cauchy case on the complex plane (cf. 
[32h using a partition of unity. To see (ii) note that the integral over LOp ooj^ is 
a (multivalued) holomorphic function, while the integral over the rest is a linear 

combination of the (multivalued) meromorphic functions ly respectively Ij . By 

construction, Ij has at most simple poles at the points fi and thus (ii) follows. 
Finally, to see (iii) observe that LOp 00+ restricted to S converges uniformly to zero 

as p — > 00+. Moreover, Ij (00+) = and hence (iii) holds. □ 

Now, let the operator Cw ■ -^^(S) — !• i^(S) be defined by 
(A.8) C^f^C+{fw-) + C-{fw+) 

for a 2 X 2 matrix valued /, where 

?«+ = 6+ — I and ui_ — I — 6_ . 

Theorem A. 5. Assume Huvothesis \A.l\ and let mo G be given. 
Assume that ji solves the singular integral equation 

(A. 9) = mo + Cui/i in L^iY,). 

Then m be defined by the integral formula 

(A.IO) m = mo + C{nw) on M \ S, 

where w = 10+ + W-, is a solution of the meromorphic Riemann-Hilbert problem 
Conversely, if m is a solution of (jA.l[) . then 11 defined via ^ — m±bj^ solves 

Proof. Suppose ji solves (|A.9p . To show that m defined above solves (jA.ip note 
that 

m± —1 + C± (nw) . 
Thus, using C+ — C_ = I and the definition of C^, we obtain 
m+ = (mo + C+{fiw)) = (mo + C+{^w+) + C+(/xw_)) 

~ (mo + + C-{fiw+) + C+{^W-)) = (mo + + C^^) 

= + 

and similarly m„ = /i(I — w_). Hence m+b^^ — fj, — ni-bz} and thus m+ — 
m_ (&_ ) ^ ^ 6+ . This proves the jump condition. That m has the right devisor and the 
correct normalization at 00+ follows from Theorem lA. 41 fii) and (iii), respectively. 

Conversely, if m is a solution of the Riemann-Hilbert problem (jA.l[) . then we can 
set /i — m+b^^ = m_6l^ and define m by (jA.lOp . To see that in fact m — rh holds, 
observe that both satisfy the same additive jump condition m+ — m_ = 7714. — m_ = 
fj,w. Hence the difference m — rh has no jump and thus must be meromorphic. 
Moreover, by the divisor conditions (mi — mi) > —Vf,* and (m2 — m2) > —Vp,, 
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the Riemann-Roch theorem imphes that m — m is constant. By our normahzation 
at 00+ this constant must be the zero vector. Thus m — fh and as before one 
computes 

m+ = - /i + niQ + CnifJ., 
showing that (jA.9p holds. □ 

Remark A. 6. (i). The theorem stated above does not address uniqueness. This 
will be done in Theorem \B.1\ under an additional symmetry assumption. 

(ii). The notation b+,b- is meant to make one think of the example = 
{b^)^^b^ in Section^ but the theorem above is fairly general. In particular it also 
applies to the trivial factorizations = IJ"^ = J^l. 

We are interested in the formula (jA.10|) evaluated at C)o_. We write it as 
m(oo_) = (too + C{^w)){oo-) 

(A.ll) 



Too 



{I-Cw) ^(TOo)wfl^_ 



and we perturb it with respect to w while keeping the contour E fixed. 

Hence we have a formula for the solution of our Riemann-Hilbert problem m[z) 
in terms of (I— C„,)~^too and this clearly raises the question of bounded invertibility 
of I — Cw. This follows from Fredholm theory (cf. e.g. 

Lemma A. 7. Assume Hypothesis \A.1[ Then the operator I — Cw is Fredholm of 
index zero, 

(A.12) ind(I- Cu,) = 0. 

Proof. Using the Bishop-Kodama theorem |26j we can approximate w± by func- 
tions which are analytic in a neighborhood of E and hence, since the norm limits 
of compact operators are compact, we can assume that w± are analytic in a neigh- 
borhood of E without loss of generality. 
First of all one can easily check that 

(A.13) (I - C„)(I - C_„) = (I - -C^)=I- T„, 

where T^if) = C-[C-{fw+)w+\. But Tu,(f) is a compact operator. Indeed, sup- 
pose /„ g L'^C^) converges weakly to zero. We will show that ||Tu,/„||2,2 — 0. 

Using the analyticity of w+ in a neighborhood of E and the definition of C_ , we 
can slightly deform the contour E to some contour E' close to E, on the right, and 
have, by Cauchy's theorem, 

(A.14) r^/„(p) = ^ / {CUnW+)w+)^i. 

Now clearly (C(/„u)+)ui+)(p) as n — > 00. and since also |(C(/„w+)'u;+)(p)| < 
const ||/„||i2||i/;+||ioo < const we infer \\Tjjj fn\\ 1^2 — )• by virtue of the dominated 
convergence theorem. 

Hence by [SH Thm. 1.4.3] I — is Fredholm. Moreover, consider ind(I — eCw) 
for < £ < 1 and recall that ind(I — eC^) is continuous with respect to e ( ^34j 
Thm. 1.3.8]). Since it is an integer, it has to be constant, that is, ind(I — C^) = 
ind(I) =0. □ 

By the Fredholm alternative, it follows that to show the bounded invertibility of 
I — Cu. we only need to show that ker(I — Cw) = 0. The latter being equivalent to 
unique solvability of the corresponding vanishing Riemann-Hilbert problem. 
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Corollary A. 8. Assume Hypothesis \A.l[ 

A unique solution of the Riemann-Hilbert problem (jA.l|) exists if and only if the 
corresponding vanishing Riemann-Hilbert problem, where the normalization condi- 
tion is given by m(oo+) = (O O) , has at most one solution. 

We are interested in comparing two Riemann-Hilbert problems associated with 
respective jumps wq and w with \\w — wo||oc small, where 

(A.15) ll^lloo = ||w+||i=o(s) + ||w_||ioo(s). 

For such a situation we have the following result: 

Theorem A. 9. Assume that for some data vUq the operator 

(A.16) I-C^* :L2(s)^L2(E) 

has a bounded inverse, where the bound is independent oft. 
Furthermore, assume w* satisfies 

(A.17) \\w* - w*o\U < ait) 

for some function a{t) — )■ as i — )■ oo. Then (I — C^^t)^^ : L'^{T,) — > -/^^(S) also 
exists for sufficiently large t and the associated solutions of the Riemann-Hilbert 
problems (jA.ip only differ by 0(Q:(t)). 

Proof Follows as in [17] • □ 
Appendix B. A uniqueness theorem for factorization problems on a 

RiEMANN SURFACE 

In the case where the underlying spectral curve is the complex plane it is often 
useful to have a theorem guaranteeing existence of a solution of a Riemann-Hilbert 
problem under some symmetry conditions. One such is, for example, the Schwarz 
reflection theorem provided in |46) . In this section we state and prove an analogous 
theorem where the underlying spectral curve is our hyperelliptic curve with real 
branch cuts. 

For any matrix (or vector) M we denote its adjoint (transpose of complex con- 
jugate) as M*. Then we have 

Theorem B.l. Assume in addition to Hvvothesis lA . 1\ assume that fij G [i?2j-i, E2j\ 
and that E is symmetric under sheet exchange plus conjugation = S J such that 

(i) Jip*) - Jip)*, forpeJ:\7r-^{a{Hg)), 

(ii) Re(J(p)) = 5(J(p) + Jip)*) is positive definite for p E n^^iaiHq)), 
(Hi) J is analytic in a neighborhood ofE. 

Then the vector Riemann-Hilbert problem (jA.ip on M has always a unique so- 
lution. 

Note here that the +- side of the contour is mapped to the — side under sheet 
exchange. In particular, the theorem holds if J = I, that is there is no jump, on 
n-\aiH,)). 

Proof. By Corollary lA. 81 it suffices to show that the corresponding vanishing prob- 
lem has only the trivial solution. 

Our strategy is to apply Cauchy's integral theorem to 

TO(p)m*(p*) — mi(p)mi(p*) -|- TO2(p)m2(p*). 
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To this end we will multiply it by a meromorpliic differential dQ which has zeros 
at /i and fi* and simple poles at oo-t such that the differential m{p)m* (jj*)dfl{p) is 
holomorphic away from the contour. 
Indeed let 

(B.f) dn = -I ' dT: 



R 



1/2 
2g+2 



and note that ^^^^^ is a Herglotz-Nevanlinna function. That is, it has positive 

imaginary part in the upper half-plane (and it is purely imaginary on a{Hq)). Hence 
m(j))m'^ {p)dfl{p) will be positive on Tr^^{a{Hq)). 
Consider then the integral 

(B.2) / m{p)m*{p*)dn{p), 

Jd 

where D is a ^-invariant contour consisting of one small loop in every connected 
component of M\I]. Clearly the above integral is zero by Cauchy's residue theorem. 
We will deform Z? to a *- invariant contour consisting of two parts, one, say -D+, 
wrapping around the part of S lying on 11+ and the + side of TT~^{a{Hq)) and the 
other being £)_ = 1?+*. 

For each component Ej of E \ TT^^{a{Hq)) there are two contributions to the 
integral on the deformed contour: 



m+(j))m*_{p*)dfl — / m-{p)J(p)m*_{p*)dVL and 

I m^{p)ml{p*)dn = [ m_{p)J*{p*)m*_^{p*)dn. 

Because of condition (i) the two integrals cancel each other. 
In view of the above and using Cauchy's theorem, one gets 



D 



0= / m{p)m*{p*)dn 

[m+(p)m1(p*) + m_(p)TO^(p*)]drj 



7r-i(^(H,)) 

m_(p)(J(p) + J*{p*))m-ip*)dn. 

7r-i(^(H,)) 

By condition (ii) it now follows that m_ = and hence m — C{fiw) with fi = m_ = 
by Theorem lA.51 (where we used the trivial factorization 6_ = I and 6+ = J). □ 

Remark B.2. The same proof also shows uniqueness for the following symmetric 
vector Riemann-Hilbert problem on M 

m+{p) = m^(p)J{p), _p e E, 
(B.3) m{p*)^m{p)(^^ 

m(oo+) = (l *) , (mi)>-I?/i«, (m2) > -P/i 
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where J{z), S, and Vfi satisfy the same assumptions as in the previous theo- 
rem. Just note that in this case the symmetry assumption implies m{p)m*i^*) = 
mi{p)m.2{p) + m.2{p)mi{p). 
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